Properties of Ca sparks revealed by four-dimensional confocal 
imaging of cardiac muscle 
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Parameters (amplitude, width, kinetics) of Ca 2+ sparks imaged confocally are affected by errors when the spark 
source is not in focus. To identify sparks that were in focus, we used fast scanning (LSM 5 LIVE; Carl Zeiss) com- 
bined with fast piezoelectric focusing to acquire x—y images in three planes at 1-jim separation (x-y-z-t mode) . 
In 3,000 x—y scans in each of 34 membrane-permeabilized cat atrial cardiomyocytes, 6,906 sparks were detected. 767 
sparks were in focus. They had greater amplitude, but their spatial width and rise time were similar compared with 
all sparks recorded. Their distribution of amplitudes had a mode at AF/F 0 = 0.7. The Ca 2+ release current under- 
lying in-focus sparks was 1 1 pA, requiring 20 to 30 open channels, a number at the high end of earlier estimates. 
Spark frequency was greater than in earlier imaging studies of permeabilized ventricular cells, suggesting a greater 
susceptibility to excitation, which could have functional relevance for atrial cells. Ca 2+ release flux peaked earlier 
than the time of peak fluorescence and then decayed, consistent with significant sarcoplasmic reticulum (SR) de- 
pletion. The evolution of fluorescence and release flux were strikingly similar for in-focus sparks of different rise 
time ( 7) . Spark termination involves both depletion of Ca 2+ in the SR and channel closure, which may be synchro- 
nized by depletion. The observation of similar flux in sparks of different T requires either that channel closure and 
other termination processes be independent of the determinants of flux (including [Ca 2+ ]sii) or that different 
channel clusters respond to [Ca 2+ ] SR with different sensitivity. 



INTRODUCTION 

In striated muscles, action potentials cause intracellular 
Ca 2+ channels to open, and the ensuing Ca 2+ release ini- 
tiates contraction. In cardiac muscle and in skeletal 
muscle of some taxonomic classes, Ca 2+ release is com- 
posed of discrete events from SR Ca 2+ channels (RyRs) 
of fairly stereotypical appearance, which have been visu- 
alized with confocal fluorescence microscopy and fluo- 
rescent Ca 2+ -sensitive dyes, and termed "Ca 2+ sparks" 
(Cheng et al., 1993; Nelson et al., 1995; Tsugorka et al., 
1995). Ca 2+ sparks are a collective phenomenon pro- 
duced by a group of mutually interacting RyR channels 
within a couplon (Stern et al., 1997, 1999), defined as 
the release channels in a cluster or Ca 2+ release unit 
(Franzini- Armstrong and Jorgensen, 1994) together 
with its associated L-type Ca 2+ channels and/or voltage 
sensors. Ca 2+ sparks occur in ventricular (Niggli and 
Shirokova, 2007) as well as in atrial myocytes (e.g., Blatter 
et al., 1997; Kockskamper et al., 2001; Sheehan et al., 
2006) . In ventricular myocytes, the well-developed 3-D 
network of transverse (t) tubules (Soeller and Cannell, 
1999) ensures that in response to an action potential, 
all couplons are activated synchronously, which results 
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in a highly uniform Ca release throughout the entire 
cell volume (Cheng et al., 1994). In contrast, in atrial 
myocytes, the t-tubular network is poorly developed or 
entirely absent (Htiser et al., 1996; Cordeiro et al., 2001; 
Mackenzie et al., 2001), and close apposition of surface 
membrane Ca 2+ channels and RyRs only exists in the 
cell periphery. The SR, however, extends throughout 
the entire cell and contains RyRs capable of Ca 2+ release 
and generation of Ca 2+ sparks (Sheehan et al., 2006). 
Although sparks appear to require the phenomenon of 
CICR (Endo et al., 1970; Fabiato and Fabiato, 1978; 
Cheng et al., 1993; Klein et al., 1996) for both initiation 
and propagation within the channel cluster, many other 
mechanistic details of these events remain unknown or 
poorly understood. Thus, numbers of channels involved, 
time course of their unitary current, duration of open- 
ing of the individual channels within a couplon, the 
mechanisms underlying their closure, the degree of local 
intra-SR depletion, and even the behavior of channels 
and currents at the time when the spark peaks have 
remained uncertain or are just starting to be elucidated. 

Many of these unknowns could be eliminated by an 
accurate quantification of sparks, in terms of their so- 
called morphometric parameters. These include peak 
amplitude, spatial width (or full width at half maximum 
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[FWHM]), rise time, decay time constant and other 
measures of duration, and integral evaluations of inten- 
sity or signal mass (ZhuGe et al., 2000; Chandler et al., 
2003) . For instance, the local amplitude of the change 
in [Ca 2+ ] together with knowledge of the unitary chan- 
nel current will provide an indication of the number of 
channels open. Similar insights will be derived from the 
accurate knowledge of spatial width. In this regard, a 
combination of modeling and experimental testing 
have shown that spatial width of sparks grows with the 
time that the release channels remain open (Zhou et al., 
2003), but widths that reach values beyond certain 
limits at early times, for example in "protoplatykurtic" 
sparks (Zhou et al., 2005) , reflect spatially large sources 
rather than small groups of channels. Finally, detailed 
knowledge of spatial aspects of sparks and their evolu- 
tion in time can be used to reconstruct the underlying 
release current by different methods that yield mutually 
consistent results (see Rfos et al., 1999, and Soeller and 
Cannell, 2002, for alternative methods, named "back- 
ward" and "forward" to distinguish whether the effect 
[the change in [Ca 2+ ]] or the cause [the release current] 
is the starting point of the calculation). An accurate 
morphometry of sparks is therefore essential for the 
elucidation of their control mechanisms. 

Crucial determinants of control of Ca 2+ sparks are 
believed to depend on cytosolic Ca 2+ itself, including 
complementary mechanisms of activation (CICR) and 
inactivation (Ca 2+ -dependent inactivation; Sham et al., 
1998, for cardiac; Baylor et al., 1983, and Melzer et al., 
1984, for skeletal). These result in channels that do 
not open or close independendy but gate concertedly. 
Whether or not channels gate independendy changes 
another morphometric property, the distribution of 
spark amplitudes, which will be a sum of decaying expo- 
nential functions of amplitude if arising from openings 
of individual channels that are mutually independent 
(and are "Markovian" or memoryless, e.g., Colquhoun 
and Hawkes, 1983; see also http:/ / en.wikipedia.org/ 
wiki/Markov_process for a brief definition) , but will typi- 
cally have a mode or preferred amplitude if channels 
gate concertedly (Bridge et al., 1999; Cannell and 
Soeller, 1999; Rfos et al., 2001; Wang et al., 2002). For 
instance, it has been argued that the rise time character- 
istics are incompatible with the idea that a spark arises 
from a single RyR channel with a reversible Markovian 
gating scheme (Shirokova et al., 1999; Rfos et al., 2001; 
Wang et al., 2002). On the other hand the spark ampli- 
tude experimentally recorded in line-scan mode always 
obeys a monotonically decaying distribution, regardless 
of the true spark amplitude (Cheng et al., 1999; Rfos 
et al., 2001). 

A major obstacle to reaching the desired accuracy in 
measures of spark parameters is that measurements are 
done on confocal images, including line scans (x-t, 
herein called 2-D) and x-y scans (referred to as 3-D 



imaging). In either case, sparks are imaged without 
knowledge of the location of its originating couplon 
relative to the line or plane of scanning. Starting with 
the work of Pratusevich and Balke ( 1996) , it has become 
increasingly clear that all morphometric spark parame- 
ters are affected by the out-of-focus error. Out-of-focus 
errors include reduced amplitude, increased spatial 
spreading, and slower kinetics. Among the parameters, 
amplitude is most affected, decaying sharply with distance 
from source to scanning line or plane, whereas rise time 
and FWHM are less distorted (Pratusevich and Balke, 
1996; Smith et al., 1998) . It is because of this susceptibility 
to focus error that the distribution of amplitudes of sparks 
detected in line scans must be monotonously decreasing. 

Theoretical approaches have been developed to recover 
the true spark amplitude, allowing quantitative correc- 
tions of spark amplitudes as if all sparks would have been 
in focus (Shirokova and Rfos, 1997; Izu et al., 1998; Rfos 
et al., 2001). Specifically Rfos et al. (2001) provide an 
equation to derive the distribution of true amplitudes 
from the distribution histogram of measured amplitudes. 

Although these procedures have value, as they de- 
monstrated, for example, a mode in the distribution 
of true amplitudes, consistent with expectations for 
groups of nonindependent channels, the correction 
formula has not been verified. Moreover, it requires 
knowledge of other parameters, which are also affected by 
the out-of-focus error and therefore are not knowable 
with certainty. 

Here, we take advantage of two novel tools, the fast x-y 
scanning provided by a slit confocal scanner (LSM 5 LIVE; 
Carl Zeiss) and the ability to move the plane of focus verti- 
cally in rapid and reproducible manner to achieve x-y-z-t 
(or 4-D) imaging of sparks. 4-D scanning eliminates the 
out-of-focus error because it allows identification of 
sparks that are imaged in focus. In this way, we charac- 
terize morphometrically a large group of sparks in 
focus. We use these, for the first time well-determined 
measures, to both clarify aspects of control of Ca 2+ re- 
lease and produce a strict test of the theory of scanning 
in lower dimensions. 

Part of this work has been presented in abstract form 
(Shkryl, V.M., and LA. Blatter. 2008. Biophysical Society 
52nd Annual Meeting. Abstr. 495; Shkryl, V.M., and LA. 
Blatter. 2009. Biophysical Society 53rd Annual Meeting. 
Abstr. 1415; Shkryl et al. 2011. Biophysical Society 55th 
Annual Meeting. Abstr. 3032) . 

MATERIALS AND METHODS 

Cell isolation and solutions 

Single myocytes were isolated from cat atria as described previ- 
ously (Wu et al., 1991; Kockskamper and Blatter, 2002; Sheehan 
and Blatter, 2003). The procedure for cell isolation was fully 
approved by the Institutional Animal Care and Use Committee. 
In brief, adult mongrel cats of either sex were anesthetized with 
sodium pentobarbital (50 mg kg -1 ). After thoracotomy, hearts 
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were quickly excised, mounted on a Langendorff apparatus, and 
retrogradely perfused with oxygenated collagenase-containing 
solution at 37°C. Myocytes were used 1-6 h after isolation. Freshly 
isolated myocytes were plated on glass coverslips in normal 
Tyrode's solution (composition in mM: 135 NaCl, 4 KC1, 2 CaCl 2 , 
1 MgCl 2 , 10 D-glucose, and 10 HEPES, with pH adjusted to 7.4 
withNaOH). 

Imaging of Ca 2+ sparks was performed in cells permeabilized 
with saponin (0.005%) for 40 s (Zima et al., 2003) in an internal 
solution composed of (in mM): 120 potassium aspartate, 15 KC1, 
5 KH 2 P0 4 , 5 MgATP, 0.35 EGTA, 0.14 CaCl 2 , 0.75 MgCl,, 10 phos- 
phocreatine sodium salt, 4% dextran, and 10 HEPES, with 
pH adjusted to 7.2 with KOH. Subsequendy, the bath solution was 
exchanged to a saponin-free internal solution containing 40 uM 
fluo-4 pentapotassium salt (Invitrogen) . The free [Ca 2+ ] and [Mg 2+ ] 
of this solution were 100 nM and 1 mM, respectively. All experi- 
ments were performed at room temperature (22-24°C). Chemi- 
cals were purchased from Sigma-Aldrich. 

Fluo-4 was excited at 488 nm, and emission was collected at 
>520 nm. The cell was positioned with its long axis parallel to the 
x scanning direction (which in slit scanners is acquired simultane- 
ously at all points of a 512-pixel line). 

4-D confocal image acquisition and multidimensional analysis 
of sparks 

Fluorescence image acquisition was performed on a slit scanning 
confocal microscope (LSM 5 LIVE; Carl Zeiss) equipped with a 
63x, 1.20 n.a. water-immersion objective (C-Apochromat; Carl 
Zeiss) and a piezoelectric focusing attachment, which by moving 
the objective allowed us to collect three x-y images (512 x 31 pix- 
els; pixel distance of 0.42 um) in succession, separated by 1 um in 
the vertical (z) direction. The point-spread function of the system 
with the focusing attachment had FWHM values of 0.52, 0.46, and 
1.25 pm, respectively, in the x, y, and z directions. The scan time 
for an individual image was 1.16 ms (line dwell time of 15.7 ps), 
and the acquisition time for one trio of images was 5.6 ms (Fig. 1 A; 
note that the acquisition time for a trio is longer than the sum 
of scan times for three images because of the time required for 
scanner fly back and movement of the piezoelectric focusing 
attachment). A total of 1,000 trios, or 3,000 x-y images, were ac- 
quired for each cell at the same location. To refer to such series of 



images we will use interchangeably the terms "x-y-z-t series" and 
"4-D series." We will refer to series at a single z position as "x-y-f 
or "3-D series." Similarly, "x-t series" will be also called "2-D" or 
line scans. 

Images in 4-D and 3-D series were normalized to resting fluo- 
rescence by a procedure in which individual images F(x,y) are di- 
vided pixel by pixel by a normalizing image F 0 (x,y) constructed by 
averaging every image with spark regions blanked. The set of 
blanked images was also vised to compute a bleaching correction 
function (see Brum et al., 2000, for details). 

Sparks were detected automatically in normalized, bleach- 
corrected x-y-t series using the algorithm developed by Brum et al. 
(2000), which is an adaptation to 3-D scans of the procedure in- 
troduced for line scans (Cheng et al., 1999). The procedure iden- 
tifies the location coordinates of the center of mass of each spark, 
a maximum of &F/F 0 (where \F= F — F 0 ) in each image, a sub- 
series of consecutive images recognized as corresponding to the 
same spark and a peak of the maximum of AF/F 0 within the sub- 
series, which we refer to as peak amplitude or just amplitude (a). 
The analysis also produced for each spark a time of first detec- 
tion, an interval between first detection and the attainment of 
peak amplitude (rise time, T), plus several measures of the spatial 
properties, including width in x and y, and eccentricity. 

For this study, we specifically developed a definition and auto- 
matic determination of FWHM. We defined FWHM as the diam- 
eter of a circle having the same area as the 2-D intercept of the 
spark at the time of peak amplitude and the plane located at half 
this amplitude (i.e., the x-y region where AF/F 0 was less than half- 
amplitude) . At every point in time during its evolution, the spark 
fluorescence is a function of x and y that resembles a bivariate 
Gaussian (examples in Fig. 1 A). An alternative set of parameters 
was derived by fitting to every spark in every image a bivariate 
Gaussian function. Parameters of fit included two coordinates of 
the center, an amplitude, and elliptic semi-axes (a and b) . In this 
description, the FWHM is equal to the geometric mean ((2a x 
2b) 1/2 ) of the two elliptic axes. Both definitions of parameters 
yielded essentially identical results. 

Other numerical routines were used to produce subsets of the 
4-D database, including 3-D (x-y-t) series at constant zand 2-D (x-t) 
series at constant z and y. In turn, these simpler subsets were 
analyzed for event detection and flux calculation by programs 
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Figure 1 . Scanning in three dimensions of space 
and time. (A) Square sections of successive x-y scan 
images of fluorescence of an atrial cardiomyocyte. 
As illustrated in the diagram, the vertical positions 
correspond to the z positions of three scanned 
planes (zi = — 1 pm, z 2 = 0 pm, and z 3 = +1 pm). As 
indicated by the horizontal arrows, the acquisition 
time for a trio of x-y scans was 5.6 ms. (B) Profiles 
(AF(x) / F 0 (x)) at the central y value of the images in 
A, color-coded to z value. Asterisks mark image trios 
where the amplitude was largest in z 2 . The example 
spark appears to be in focus. 
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developed earlier (Cheng et al., 1999; Brum et al., 2000). All of 
the numerics were implemented in the IDL programming en- 
vironment (ITT Visual Information Solutions). 

Methods to decide whether a spark is in focus 

A spark that appears to be in focus is illustrated in Fig. 1. In Fig. 1 B, 
fluorescence intensity profiles (2 pixels wide) of intensity traced 
through the centers of the sparks, of which all correspond to the 
same location in successive images (Fig. 1 A) , are shown. It can be 
seen that in several trios of images, the peaks of fluorescence in 
the zg plane are greater (asterisks) than in the Zi and z 3 planes. 
This spark therefore satisfies the simplest criterion for being in 
focus, namely, that its amplitude be greatest in the central plane, 
Z2. This simple criterion, however, failed for several reasons. It is 
strongly biased toward acceptance of sparks of narrow width and 
is susceptible to limitations in scanning speed, limitations that will 
become clearer when describing Fig. 5. After testing several other 
procedures, we opted for a "three-point criterion," so named be- 
cause it uses asymmetrically all three measures of amplitude, from 
all three focal planes (z 1; z%, and z 3 ). The criterion, which is de- 
scribed in the Appendix, allows defining the 2-D vertical thickness 
of the slice where the sparks must originate to be considered in 
focus. All sparks deemed to be in focus in this study were identi- 
fied by this criterion, with 2D = 0.2 pm. 

Reconstruction of line-scan images from 3-D scans 

One of the aims of this work is to compare the properties of 
sparks imaged in focus with those inferred from images obtained 
in line-scan mode, which include sparks that may or may not be 
in focus. For this purpose, we derived from the database of 4-D 
series a set of images equivalent to line scans. The extraction of 
such reconstructed line scans is illustrated with Fig. 4 A. 

For each image at a single value of z (z 2 = 0 pm), a line was ex- 
tracted, namely the group of 512 values F(xj, y n , z?) at a set value of 
y (% = 6 P m > the approximate center of the range of y values), 
with j varying from 1 to 512. The lines extracted at the same y and 
z values in 1,000 successive images constitute together a 2-D series 
F( jij, i k ) , the reconstructed line-scan image (with k varying from 1 
to 1 ,000) . One such line scan is represented in Fig. 4 A. The interval 
between sets of x-y-z scans (or the interval between x-y trios) was 
5.6 ms, which is then the temporal increment of the line scan. 



RESULTS 

Cat atrial myocytes, which are devoid of a t-tubular system 
(Hiiser et al., 1996), are known to generate spontane- 
ous Ca 2+ sparks from RyR clusters located in the non- 
junctional SR (nj-SR) membranes (Sheehan et al., 2006) . 
Thus, atrial Ca 2+ sparks from nj-SR result entirely from 
SR Ca 2+ release and lack any potential component from 
Ca 2+ influx. We studied Ca 2+ spark properties in these 
cells using a fast slit confocal scanner (cf. Toomre and 
Pawley, 2006, for brief description) in combination with 
a fast piezoelectric focusing attachment to acquire con- 
focal x-y images successively in three vertical planes, at 
1-pm separation. The imaging process, which we call 
4-D scanning, is illustrated in Fig. 1 (see Table 1 for a list 
of terms and symbol definitions). Fig. 1 A includes a 
series of images from an atrial cardiomyocyte (with 
plasma membrane permeabilized and in a solution con- 
taining fluo-4). Square regions from the successive im- 
ages are shown, including a spark that appears to be in 
focus. The relative position in the figure of each square 



image in the horizontal direction suggests the time at 
which it was acquired, whereas the vertical positions in 
the figure correspond to the three z positions of the 
scanned planes {z\ = — 1 pm, z» = 0 pm, and z 3 = +1 pm) . 
Fig. 1 B plots the central profiles (A.F(x) / F 0 (x) , where x 
corresponds to the longitudinal axis of the cell) through 
every frame, color-coded to z position. 

The example spark shown in Fig. 1 appears to be in 
focus because it is somewhat more intense in the central 
(z 2 ) plane of scanning (asterisks in B mark image trios 
where the amplitude was largest in za). As stated in 
Materials and methods, the simple criterion that an in- 
focus spark would be detected with greater amplitude 
in the central plane led to biased results. Therefore, the 
three-point criterion (cf. Appendix), which is also satis- 
fied by this example, was used to identify in-focus sparks. 
In the following sections, we will explore the properties 
of sparks in focus and compare them with those of all 
detected sparks. 

Properties of in-focus sparks 

An algorithm of automatic detection developed for x-y 
images (Brum et al., 2000) was applied to 34,000 x-y-z 
scans (i.e., trios of x-y images at three z values) obtained 
from 34 cells at 1,000 x-y-z (or 3,000 x-y) scans per cell. 
In these, 6,906 sparks were detected, of which 767 were 
found to be in focus by the three-point criterion. 

The histograms of amplitudes (defined as a = peak 
AF/F 0 ) of the two groups are represented in Fig. 2 A. 
Note that the histogram of all sparks detected (Fig. 2 A, 

TABLE 1 

List of symbols 
Symbol Description 



x, y, z, 
*i° 

F F„, F/F„, AF/F„ 
a, FWHM, T 
g(a), 1(a), u(a) 

C,(a), L(a), U(a) 

h(z,<r) 

I 

8 

C3 

1) 
P 

[C; 
N, P„, i 



Spatial coordinates {x, horizontal, axial; y, 

transverse; z, vertical) 
Polar spatial coordinates {radial and angle) 
Fluorescence, initial, normalized, increase 
Spark amplitude, FWHM, rise time 
Probability density of the distributions of 

amplitude (in focus, 3-D, 2-D) 
Histograms, sampling of the distributions g{a), 

l{d), u{a) 

Gaussian function of vertical coordinate z with 

standard error <x 
Vertical separation of spark source and focal 

plane 

Normalized difference in fluorescence {defined 
by Eq. 1) 

"Three-point criterion" {limit value of 8 for sparks 

in focus) 
Half depth of slice in focus, pm 
Pearson's linear correlation coefficient 
Free Ca 2+ concentration, in cytosol, in SR 

Number of channels, channel open probability, 
unitary current 
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dashed black line) , referred to as L(a) , is monotonically 
decreasing, except for a region of the domain at low 
values of a (<0.4) , where detection fails. In contrast, the 
histogram of in-focus spark amplitudes (Fig. 2 A, G( a) , 
solid green trace) features a mode at about a = 0.7. The 
implications of the distribution of amplitudes of sparks 
in focus will be fully considered in Discussion. 

In Fig. 2 B, we compared in the same way as in Fig. 2 A 
the spatial width (FWHM) of detected sparks. In this 
case, the differences between the values reported for 
in-focus sparks and those of all detected sparks were 
much smaller; both distributions spread asymmetrically 
around a mode, close to 2.6 pm, but the distribution of 
in-focus sparks was slightly narrower. Clearly, the mea- 
surement of spatial width is much less compromised by 
out-of-focus error than that of spark amplitude. 

In Fig. 2 C, histograms of rise time T (from the time 
of first detection to that of peak amplitude) are shown. 
Twas measured always in the central plane (22). This 
resulted in a temporal resolution of 5.6 ms per image 
and consequently a coarse distribution. It is also neces- 
sary to point out that rise time here is defined as start- 
ing from the frame of first detection. This definition 
underestimates the true rise time and will warrant a re- 
interpretation of measured times (compare below). 

These cells produced on average rather wide sparks 
of moderate amplitude. Averages of these parameters 
are listed in Table 2 for the 4-D, 3-D, and 2-D datasets. 

The newly found method to identify sparks in focus 
obviously allows for a variety of other studies. Wang 
et al. (2004) found a quantized distribution of the rates 
of rise of fluorescence, indicative of quantized levels of 
flux, in sparks of peripheral couplons located in focus 
by activation of ventricular myocytes via loose patch 
clamp. To compare, and even though the present imag- 
ing was done at a much coarser temporal and spatial 
resolution, we computed initial rates of rise of fluores- 
cence (rate of rise in a 5.6-ms interval from first detection) 
for individual sparks in the in-focus group. The values 
had a broad distribution between 0.01 and 0.12 ms -1 . 
There was no evidence of quantization. 

Morphometric parameters of in-focus sparks 
are weakly correlated 

In addition to the morphometric parameter values, 
other properties of sparks with interesting mechanistic 
implications are the correlations among these param- 
eters. The correlations found in published studies 
have been poor for sparks of both skeletal and cardiac 
muscle cells. For example, linear correlation coeffi- 
cients (Pearson's p) between rise time and amplitude 
of 0.09 and 0.12 were found in two studies of frog skel- 
etal muscle (Rfos et al., 1999, 2008). In cardiac myo- 
cytes, correlations between amplitude, rise time, and 
width revealed p values between 0.07 and 0.17 (Shen 
etal., 2004). 



A chief factor determining the poor correlations is 
the contribution to the variance in amplitude by the 
out-of-focus error (i.e., the absence of knowledge of the 
separation £ between scanned line and spark source), 
which amounts to between 80 and 90% of the total vari- 
ance (Rfos et al., 2008) . The newly found ability to iden- 
tify sparks in focus should remove the uncertainty of 
location, therefore reducing the out-of-focus contribution 
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Figure 2. Properties of sparks in focus. Histograms of ampli- 
tudes (A), FWHM (B), and rise time (C) for all sparks detected in 
x-y scans (black, dashed), sparks in focus (green), and sparks de- 
tected in reconstructed line scans (red). Note that the amplitude 
histograms differ greatly among the three groups of sparks, but 
the other distributions diverge less. In-focus sparks have greater 
amplitudes and a mode in the distribution of this variable. Their 
modal width is not very different than that of the other groups, 
and their rise times are nearly identical to those of the entire 
set of detected sparks. The distribution of amplitudes measured 
in reconstructed line scans is heavily biased toward low values, 
whereas those of widths and rise times are more disperse than the 
corresponding ones in x-y scans. Relevant parameters of these 
distributions are listed in Table 2. 
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TABLE 2 





Avstcl^s wtoTfthowtatnc j 








Mode 


Amplitude 


Rise time 


FWHM 


N 










sparks 






ms 


y/m 




4-D 


1.05 


13.5 


3.35 


767 


In-focus 


(0.01) 


(0.24) 


(0.02) 




3-D 


0.82 


10.7 


3.18 


6,906 




(0.005) 


(0.10) 


(0.01) 




2-D 


0.61 


17.4 


3.85 


2,910 




(0.007) 


(0.22) 


(0.03) 





The SEM is in parentheses. The row with heading "4-D In focus" lists values for sparks detected in x-y-z-t series, the 4-D mode, which satisfy the in-focus 
criterion, x-y-z-t series were obtained from 34 cells; each series consists of 1 ,000 trios of x-y images successively at three z positions, for a total of 3,000 images 
per series. The row with heading "3-D" includes all sparks detected in x-y-t series at the zq position (z - 0) , whether or not they were in focus. "2-D" refers to 
sparks detected in line-scan images constructed from the x-y-t series at position z<> as described in Materials and methods. Fig. 4 has an example line-scan 
image constructed this way. It should be noted that the three modes (4-D, 3-D, and 2-D) are just different ways to analyze the full 4-D database. 



to the variance in amplitude, and thus provide a better 
gauge of existing correlations. 

The correlations among the parameters amplitude, 
rise time, and FWHM are illustrated by joint histograms 
in Fig. 3. The left-side panels represent the histograms 
of all sparks (3-D mode), whereas the corresponding 
histograms of events in focus are on the right side. Positive 



correlations were found in every case and were usually 
modest. Amplitude and width had p = 0.04 for all 
sparks and 0.14 for those in focus. For rise time and 
width, p increased from 0.14 to 0.26 in going from all 
sparks to those in focus, whereas for rise time and am- 
plitude, p increased from 0.33 to 0.36. In sum, correla- 
tions between all three morphometric parameters were 
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Figure 3. Correlations between spark parameters. 
Joint histograms of amplitude and width (A and B), 
rise time and width (C and D), and rise time and 
amplitude (E and F) for the set of all sparks detected 
in x-y images (A, C, and E) and for in-focus sparks 
(B, D, and F). The total number of sparks included 
was 6,906 (all sparks) and 767 (in-focus sparks). The 
first-order correlation coefficient p (a, FWHM) was 
0.037 for all sparks and 0.14 for sparks in focus, 
p (T, FWHM) was 0.14 (all) and 0.26 (in focus), 
p (T, a) was 0.33 (all) and 0.36 (in focus). The color 
table starts at 0 (black) for all panels, whereas the 
maximum (yellow) corresponds to 143, 35, 201, 85, 
149, and 70 events, respectively, for A-F. 
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positive and modest, and all of them increased for sparks 
in focus. The strongest positive correlation was found 
between rise time and amplitude. The issue of correla- 
tions will be encountered again when considering prop- 
erties of averages of sparks. 

The most relevant advantage of identifying sparks in 
focus is that the local Ca 2+ transient, determined with 
minimal error, can then be used to derive the underly- 
ing flux of Ca 2+ release. This will be done in the Discus- 
sion section. 

The morphometric properties of sparks in line-scan images 
It seemed interesting to compare the presumably true 
properties of sparks revealed by the present technique 
with those found using the line-scan (or 2-D) mode, 
which is still the standard of spark acquisition. A set of 
images equivalent to line scans was derived from the 
database of x-y-z-t series as described in Materials and 
methods and analyzed to derive their morphometric 
parameters in the standard way. One reconstructed line 
scan is shown in Fig. 4 A. A total of 2,910 sparks in 2-D 
were derived from the 4-D database. The number is less 
than in the x-y-t or 3-D case (6,906) because the method 
used one line per x—y image at a fixed z value and there- 
fore missed a substantial portion of the images, result- 
ing in fewer identified sparks. From these numbers, and 
considering that the length of scanning inside cells was 
on average 117 um, a frequency of events can be calcu- 
lated as /= 2,910 x (117pm) _1 x (5.6 s)" 1 x (34 cells)" 1 x 
100= 13.1 events xs" 1 x (100 pm)" 1 . 

The histograms of amplitude, rise time, and spatial 
width of this group are represented in red traces in Fig. 2. 
The histogram of amplitude (referred to as U(d)) decays 



monotonically (except in a region of low amplitude, 
a < 0.4, where density increases as a result of increased 
detectability) and with positive curvature, in agreement 
with expectations from the theory of line-scan imaging 
(Rios et al., 2001, and Discussion). The histogram of 
FWHM is broader for the 2-D case (Fig. 2 B, red), but 
other properties (like mode and skewness) are robusdy 
conserved in all dimensional approaches (i.e., 2-, 3-, and 
4-D). A similar broadening occurs with the histogram of 
rise times (Fig. 2 C). 

DISCUSSION 

This work evaluates quantitatively the properties of 
sparks imaged in focus, based on the availability of a fast 
confocal scanner, which combined with a rapid change- 
of-focus mechanism allows for actual scanning in three 
spatial dimensions. This mode of scanning can be re- 
peated at a rate sufficient to characterize the evolution 
of sparks in space and time. 

4-D scanning can be used to define sparks in focus 

The identification of in-focus sparks reliably and repro- 
ducibly hinges critically on a robust criterion for the 
identification of the optical plane from where the sparks 
originate. Several approaches have been proposed. Sparks 
elicited by sparklets, the highly localized Ca 2+ signal origi- 
nating from Ca 2+ entry through voltage-gated Ca 2+ chan- 
nels in the cell membrane near a patch electrode, are 
considered to be in focus (Wang et al., 2001, 2004), 
whereas the simultaneous recording of cytosolic sparks 
and the corresponding intra-SR Ca 2+ depletion signal 
(termed Ca 2+ blinks; Brochet et al., 2005) is an objective 
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Figure 4. Properties of in-focus sparks derived from 
line-scan images. (A) Line-scan image reconstructed 
from 1 ,000 trios of x-y images obtained successively, 
at 5.6-ms intervals, from the same cell. The image 
is formed by one line of pixels taken from the cen- 
tral (%) image of every trio, at 512 values of x and 
one value of y (y = 6 um) . Size of image was reduced 
(from 512 x 1,000 to 298 x 1,000) by eliminating 
pixels corresponding to x values outside the cell 
boundary. Lines are placed vertically in the figure 
and stacked from left to right. One reconstructed 
line scan was derived for each of 34 cells studied. 
(B) Red solid trace, U(a): histogram of amplitudes 
of 2,910 sparks detected by the conventional algo- 
rithm from all .34 reconstructed line scans. This plot 
is also in Fig. 2 A. Dashed trace, g(a): distribution of 
"true" amplitudes derived substituting U(a) for u(a) 
in Eq. 1. Green trace, G(a): amplitude histogram of 
sparks in focus, reproduced from Fig. 2 A. 
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criterion that circumscribes the location to the smaller 
volume where the blink signal can be detected (Zima 
et al., 2008a,b, 2010) . Here, we demonstrate a novel ap- 
proach, which relies on the near-simultaneous imaging 
of Ca 2+ sparks in multiple planes and a robust criterion 
to determine where the spark originates. 

In principle, a straightforward criterion for deciding 
whether sparks are in focus would be a fluorescence sig- 
nal greatest in the central (22) plane. This criterion was 
not useful for several reasons. It was biased toward the 
detection of narrow sparks, and it did not take into con- 
sideration the fact that spark amplitude changes during 
the time it takes to acquire one x-y-z stack. Another limi- 
tation is that the criterion, consisting in simply comparing 
pairwise the amplitudes at z% and the other two planes, 
does not use all the information available. An alternative 
approach was based on the realization that separation 
from the focal plane, in either direction, resulted in asym- 
metric and predictable changes in the measurements 
at Zi, Zg, and Z3. In other words, if the center of a spark 
moves away from the central plane, the measured am- 
plitude at the central plane decreases, but the other two 
measures, at Z\ and Z3, change in opposite directions. 
These simple considerations, together with the fact that 
sparks are extremely well described by Gaussian func- 
tions of space, were combined in the three-point criterion 
(Fig. Al in Appendix) . The criterion had the additional 
advantage of an adjustable stringency associated with a 
measure of the depth of the volume where events would 
be identified as being in focus. 

In principle, it is not possible to know whether or not 
a selection process works correctly. This is simply be- 
cause there is no known comparable population of 
in-focus sparks. We used instead two properties expected 
of a good selection. One is that the quantitative proper- 
ties of the group of sparks classified as in focus should 
converge as the criterion is applied with increasing 
stringency (Appendix, Fig. A2) . Additionally, the limit- 
ing properties to which sparks converge should be consis- 
tent with predictions made by the theory of line scanning 
(Shirokova and Rfos, 1997; Izu et al., 1998; Cheng et al., 
1999; Rfos et al., 2001). Reliance on this theory seems 
justified, as some of its predictions have been quanti- 
tatively confirmed with simulations (Rfos et al., 2008), 
whereas others are consistent with sparks observed under 
special conditions, which assure their being in focus 
(Bridge et al., 1999; Soeller and Cannell, 2002, 2004; 
Wang etal., 2002, 2004). 

As illustrated in Fig. A2, we found convergence of the 
histogram of amplitudes as stringency increased (within 
the limitations of the test) . Specifically, the distribution 
of amplitudes lost density at low amplitudes and increased 
it at higher amplitudes, as predicted by the theory. 
There were no major changes in rise time and spatial 
width, other than a slight narrowing of the distribution 
of spatial widths, again consistent with expectations for 



sparks in focus. In increasingly stringent applications of 
the selection procedure, we also found convergence of 
width and rise time. In the next subsection, we evaluate 
quantitatively the agreement between observations and 
predictions by the theory of line scanning. 

In-focus sparks have a preferred amplitude 

A meaningful characteristic of sparks imaged in focus 
is a histogram of amplitudes G(a) with a modal value 
(Fig. 2 A) . This is in sharp contrast with the distribution 
of amplitudes reported by line or 2-D scanning. First, 
note here the distinction between U(a), the histogram 
of amplitudes reported by line scanning, and u(a), the 
distribution that such a histogram samples. U{a) could 
be found in the present work, for the same group of im- 
ages, as the histogram of amplitudes of the sparks found 
by applying our automatic detector to line scans recon- 
structed from sequences of x-y frames at Z2. The recon- 
structed line-scan images are illustrated in Fig. 4 A, and 
the U(a) of such line scans is in Fig. 2 A (red trace). 
According to the theory of line scanning, u(a) is a 
monotonically decaying function. U(a) is consistent 
with the expectation, except for a rise explained by 
missed events at very low values of a. 
Rros et al. (2001) provide an equation, 



to derive the distribution of true spark amplitudes g(a) 
from u( a) . In Fig. 4 B, we reproduce U( a) from Fig. 2 A 
(red, solid), and plot (red, dashed) the estimate of g(a) 
obtained substituting U(a) in Eq. 1. This estimated g(a) 
compares reasonably well with G(a), reproduced from 
Fig. 2 A (green). The agreement serves both as addi- 
tional confirmation of the new method to define sparks 
in focus and as confirmation of significant properties of 
Ca 2+ sparks. 

Specifically, we confirm earlier conclusions, reached 
for both skeletal and cardiac ventricular muscle, that 
sparks have a preferential amplitude. In both tissues, the 
conclusions applied to special conditions. In the case of 
skeletal muscle, a modal amplitude was found in g(d) 
functions derived via Eq. 1 in cells exposed to low con- 
centrations of caffeine (Rros etal., 2001). In cardiac mus- 
cle, a modal G(d) was found in sparks elicited by sparklets 
near a patch electrode (therefore considered to be in 
focus; Wang et al., 2001, 2004). A modal amplitude was 
also found for groups of sparks originating from the 
same location in a cell, presumably from a single cou- 
plon firing repetitively (Bridge et al., 1999; Soeller and 
Cannell, 2002). Furthermore, the amplitudes of local SR 
Ca 2+ depletion signals corresponding to Ca 2+ sparks (Ca 2+ 
blinks) reveal a modal distribution (Zima et al., 2008b). 

As stated in earlier work, the main implication of the 
existence of a preferential amplitude of sparks is that 
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their sources are not single Markovian channels (Bridge 
et al., 1999; Ri'os et al., 2001; Wang et al., 2001, 2004). 
The existence of a mode in the distribution of ampli- 
tudes has been reproduced by models in which sparks 
are generated by channel clusters, interacting via CICR, 
calcium-dependent inactivation, or mechanical contact. 

Although the simplicity of the procedure to find sparks 
in focus and its initial tests and results are reassuring, 
there are also inherent limitations to the method. 
First, only three vertical (z) positions were scanned, and 
the spatial resolution was set low in every dimension 
(as a compromise between the increased volume of 
data and the finite rate of data acquisition) . Second, the 
vertical scanning reduced the time resolution of the re- 
cordings (in our case to 5.6 ms between successive x-y-z 
scans or image trios) , which particularly affects the 
accuracy of the recording of rapidly changing features 
of the fluorescence signal, such as the rise time. 

The properties of sparks in focus can be determined 
with precision 

Fig. 5 shows averages of 150 sparks that were found to 
be in focus and had the additional common property of 
a T= 11.2 ms. What this means is that the peak ampli- 
tude was reached two frames later than the frame of 
first detection. This number, 11.2 ms, has precision lim- 
ited by the large (5.6-ms) interval. In Fig. 5, images are 
arranged in three columns and nine rows. Each column 
contains averages at the same z position. Each row rep- 
resents images (image trio, cf. Materials and methods) 
obtained in the same 5.6-ms interval. The interval of 
first detection is labeled c. At this time, the spark be- 
came sufficiently different from noise to be accounted 
for by the detection algorithm. The time of first detec- 
tion (t = 0) occurs substantially later than the start of the 
event. This is clearly shown in Fig. 5, where the sparks 
that form the average started to develop in row a, while 
they were first detected 11.2 ms later in row c. The dis- 
crepancy between start of the event and time of detec- 
tion is the result of the chosen stringency of the spark 
detection algorithm, and leads to an underestimate of 
rise time (by ^11.2 ms) . 

It should be kept in mind that the images were ob- 
tained in the order Z\, z%, z$. This means that images on 
the right-most column were obtained later, resulting in 
a systematic excess in amplitude during the rising phase 
and a corresponding deficit during the decaying phase. 
These differences are especially noticeable at times when 
the spark is changing rapidly (as in rows b, c, f, and g) . 
At these times, the simple expectation that an in-focus 
spark should be measured with greatest amplitude at 
Z2 is clearly invalid. This property helps explain why a 
more elaborate criterion was needed to decide whether 
a spark is in focus. 

The morphometric parameter values of this average 
are fairly representative of the histograms in Fig. 3. The 



amplitude is 1.04, and the FWHM is 2.98 um. The aver- 
aged spark data (Fig. 5) were also used for the Ca 2+ 
release flux calculation discussed in the next subsection. 
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Figure 5. The average spark of T= 11.2 ms. Averages of images, 
at different times and z positions, of 150 in-focus sparks with T = 
11.2 ms. Each column of nine panels contains successive averages 
at the vertical (z) position listed at top. Each row has the aver- 
age images obtained at vertical positions Zy z%, and z 3 in the same 
5.6-ms interval. The interval of first detection is labeled c. The 
times of acquisition of the images at level z$ are listed in the third 
column. The origin of time (t= 0) is that of first detection, which 
occurs substantially later than the start of the event. Note that the 
spark average at Zg has greater fluorescence than those at z l and 
z 3 only for the trio obtained at the time of peak amplitude, row e. 
During the stage of rapid rise of fluorescence, the averages at z 3 
are greatest, whereas during the declining phase, it is position Zi 
that records the greatest waveforms. The left-side square panel in 
row e depicts the average spark at the time of peak amplitude, 
at Zg, in grayscale. The right-side square panel shows in white the x-y 
region where average F/F 0 is greater than its half-maximum value. 
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Flux of Ca 2+ release can be reliably estimated 

A key advantage of knowing that a spark is in focus is 
that the underlying flux of Ca 2+ release can then be 
derived with greater confidence. We calculated re- 
lease flux applying the "backwards" method (Rfos et al., 
1999) to averages of in-focus sparks of the same rise 
time (in the example of Fig. 5, T= 11.2 ms). The calcu- 
lation, as implemented in Rfos et al. (1999) assumes 
radial symmetry for the event. The assumption applies 
well to the present case. This is demonstrated graphi- 
cally for the average at T= 1 1 .2 ms by the black insets in 
row e of Fig. 5, which show the average at its peak in 
grayscale (left) and a mask (right) formed by the con- 
tour of the average in the region where its fluorescence 
is greater than half amplitude. This contour is an almost 
perfect circle, a symmetry that applies to the average 
but does not preclude asymmetries in the individual 
sources. Furthermore, as detailed in Materials and 
methods, two measures of spatial width, in orthogonal 
directions (x and y dimension) , were obtained for every 
spark, and this was done both by measuring diameters 
at half-maximal amplitude and by fitting a bivariate 
Gaussian to the spark. We found no significant differ- 
ence between the average values of the two orthogonal 
measures, which is of course expected if the average 
is symmetric. 

The method of release flux calculation requires as- 
signing values to concentrations and kinetic parameters 
of the relevant Ca 2+ -buffering molecules and removal 
systems, which was done as in our earlier work (Santiago 
et al., 2010). The Ca 2+ buffers and removal processes in- 
cluded in the calculation are listed, with their parame- 
ter values, in Table 3. The use of permeabilized cells is 
advantageous for flux calculations because the cytosolic 
composition is precisely known. The calculation did not 
attempt to correct sparks for the blurring effect of the 
imaging process; the effects of blurring and the deblur- 
ring correction on calculated flux were found to be 
minor by Rfos et al. ( 1999) , especially in cases where the 



sparks are spatially much wider than the point-spread 
function of the imaging system. 

The calculation and results are illustrated in Fig. 6. In 
Fig. 6 A is the fluorescence, averaged over the angle <p 
(as shown schematically in inset a), plotted as a func- 
tion of radial distance r and time from first detection. 
Inset b shows the same fluorescence in pseudocolor. 
That the fluorescence is averaged first over 150 events 
and then over the polar angle explains the near absence 
of noise in the result. 

The calculated release flux is in Fig. 6 B. The spatial 
size of the source, which can be appreciated in inset c, 
is very small compared with the fluorescence event. 

Quantitative aspects of the flux are illustrated in 
Fig. 7. The radial profiles of spark and calculated flux 
are compared in Fig. 7 A. It can be seen that the flux 
largely originates at a region of ^0.5 pm in diameter, 
which is roughly the spatial resolution of the micro- 
scope (a result suggesting that the source has dimen- 
sions below the limit of resolution). Fig. 7 B depicts 
temporal dependence of fluorescence (F/F 0 ), calculated 
free Ca 2+ concentration (pM) , and Ca 2+ release flux 
(mM/s) for the same spark average. The peak Ca 2+ re- 
lease flux was calculated as 53 mM/ s. Even though 
temporal aspects are not well described at this time res- 
olution, it can be seen that the peak of the flux pre- 
cedes the peak of free [Ca 2+ ], and both precede the 
peak of fluorescence. The graph includes a plot of Ca 2+ 
release current, calculated by integrating flux over the 
volume of the source. The current is approximately co- 
temporal with the flux and peaks at ^11 pA. This value 
of peak current is at the higher end of the range of esti- 
mates in cardiac cells (Blatter et al., 1997; Izu et al., 
2001; Cheng and Wang, 2002; Soeller and Cannell, 
2002; Wang et al., 2004; Santiago et al., 2010). Such 
level of current would require 20-30 channels simulta- 
neously passing currents of 0.3-0.5 pA (Kettlun et al., 
2003). Although large, this number of channels is still 
consistent with a source that remains below the limit 



TABLE 3 

Release flux calculation 



Binding site 


Concentration 




korF 


Diffusion coefficient 


Flux 






(fiM ms)" 1 


ms" 


pm 2 ms" 1 


mM/s 


ATP 


5,000 


Ca, 0.15 


30 


0.14 


25.3 






Mg, 2e-3 


0.2 






EGTA 


350 


4e-3 


2e-3 


0.07 


0.6 


Fluo-4 


40 


0.032 


0.032 


0.05 


1.4 


SERCA 


47 


0.5 


0.5 


N.A. 


0.9 


Troponin 


70 


5.7e-3 


0.011 


N.A. 


0.1 


Sarcolemma 


42 


0.1 


1.3 


N.A. 


0.1 


Free Ca 2+ * 


variable 


N.A. 


N.A. 


0.35 


24.7 



Model parameters (concentration, reaction rates, and diffusion coefficients) are given for all Ca 2+ ligands used in the calculation. The last column lists 
maximal contributions to the flux total by the binding of Ca 2+ to the corresponding ligand. The last row lists the free Ca 2+ term in the flux calculation 
{namely, the flux required to account for the local rate of change of [Ca 2+ ] and diffusion of the free ion). With minor changes, parameters are as in 
Santiago et al. (2010) and references therein. 
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of optical resolution. It is also consistent with estimates 
of the number of contributing channels obtained by 
noise analysis of spark amplitudes (Bridge et al., 1999), 
as well as the number of RyRs per couplon derived from 
RyR binding (Hayashi et al., 2009) and immunostaining 
studies (Soeller et al., 2007; Cannell and Kong, 2012). 
The number of channels estimated here is much greater 
than that derived from the observation of quantized flux 
in sparks from peripheral couplons (n < 6; Wang et al., 
2004). Cannell and Kong (2012) argue that a complex, 
multi-couplon source could reconcile the observation 
of quanta! sparks with the greater estimates of channel 
numbers derived in other works and confirmed here. 




The large number of simultaneously active channels 
in the present study is indicative of a greater than usual 
tendency of these channels to open. This tendency is 
manifested in the frequency of sparks detected in the 
reconstructed line scans (13.1 events; 100 um s _1 ), 
which is greater than that reported in studies of intact 
atrial and ventricular muscle (ranging from ^0.5 to 5 
events; 100 um -1 s~ h . McCall et al., 1996; Satoh et al., 
1997; Li et al., 2005; Sheehan et al., 2006; Copello et al., 
2007; Santiago etal., 2010), but is only somewhat higher 
than frequencies observed in permeabilized ventricular 
cells (3-11 events; 100 um -1 s _1 : Zima et al, 2003, 2004, 
2008a,b, 2010; Qin et al., 2009). Part of the increase 
in frequency may therefore be attributed to an effect 
of membrane permeabilization, which has been dem- 
onstrated directly for the same cell type (cat atrial 
myocytes; Sheehan et al., 2006), and part to an actual 
difference between atrial and ventricular cells (cf. dis- 
cussion below) . 

The calculation of flux permits an evaluation of the 
role of different buffer systems in shaping sparks. The 
calculation of Ca 2+ release flux produces an estimate of 
the removal flux associated with every process that is 
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Figure 6. Ca 2+ release flux of an averaged spark. (A) Normal- 
ized fluorescence, F/F 0 (r, t), obtained by averaging over the polar 
angle <p (as shown in inset a) the averaged spark of T= 11.2 ms 
(illustrated in Fig. 5). F/F 0 (r, t) is represented in pseudocolor in 
inset b. (B) Release flux, calculated from F/F 0 (r, t), represented 
as surface plot or in pseudocolor (inset c) . The "ATP component" 
(specifically, flux of Ca 2+ removal into ATP) and "free Ca 2+ com- 
ponent" (flux needed to account for local rate of change of free 
[Ca 2+ ] and diffusion of the free ion) of flux are separately plotted 
in insets d and e, respectively. Collectively, they constitute >90% 
of the total at the time of peak release flux (details in Table 3). 



Figure 7. Quantitative properties of release flux. (A) Radial de- 
pendence of fluorescence (and its proportional magnitude, con- 
centration of Ca 2+ -bound fluo-4; dashed trace) at the time of its 
peak compared with that of calculated release flux for the average 
of sparks of T= 11.2 ms. Although the FWHM of the fluorescence 
spark was 2.98 pm, that of the flux was 0.45 pm. (B) Time course 
of spatial maxima of fluorescence (green), calculated free [Ca 2+ ] 
(blue), release flux (red), and release current (black) measured 
in successive average images of sparks of T= 11.2 ms. Note that 
flux started to fall before the time of peak fluorescence, an occur- 
rence common for sparks of T> 11.2 ms. 
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taken into account. The processes and their maximal cal- 
culated removal fluxes for the average spark at T= 11.2 ms 
are listed in Table 3. Of the six processes listed, ATP 
binding and free Ca 2+ accumulation and diffusion con- 
tribute together >90% of the flux, and their contribu- 
tions, plotted separately in insets d and e of Fig. 6, have 
approximately equal maxima. 

In addition to stressing the importance of ATP in 
shaping fast cytosolic Ca 2+ transients, a notion first 
introduced by Baylor and Hollingworth (1998), this 
observation implies that the assumptions regarding 
specific properties of other ligands are unlikely to affect 
the results in any significant way. In support of this in- 
ference, alternative calculations with threefold changes 
in either direction in the diffusion coefficients assumed 
for EGTA and dye caused negligible changes in the re- 
sulting flux. Naturally, changes in the assumptions for 
ATP and Ca 2+ had a greater effect. The main conclusions, 
however, were upheld. For example, when the diffusion 
coefficient D ATP was reduced (increased) by a factor of 2, 
peak release current decreased by 8% (increased by 12%) . 
When D Cll was reduced (increased) twofold, the peak 
current was reduced by 4% (increased by 8%). Kinetic 
milestones, however, including the time to peak flux, 
and whether or not flux decreased before the peak fluo- 
rescence, were not changed in any case. 

Conversely, the concentration of free Mg 2+ , which in 
turn defines the availability of ATP, becomes an impor- 
tant factor in shaping fast local Ca 2+ transients. 

Release flux and control of channel open time 

The determination of spark properties and calculation 
of flux and current was done for spark averages differ- 
ing by rise time T. Similar averages as illustrated in 
Fig. 5 were obtained by binning in-focus sparks according 
to T, between 5.6 and 28 ms (in increments of 5.6 ms) . 
The five different values of Twill be represented by 
7], with j varying from 1 to 5. The properties of averages 
of sparks of different Tare illustrated in Fig. 8. Fig. 8 A 
plots the amplitudes at different times (with t = 0 de- 
fined as the time of first detection) . A notable aspect is 
that the early stages of all these averages are fairly simi- 
lar, with a stage of fast increase, of approximately equal 
rate for all averages, regardless of the time at which 
sparks reach peak amplitude. This common rate of rise 
indicates that the underlying flux of Ca 2+ release is simi- 
lar for the five groups (more implications of this obser- 
vation are discussed below) . 

The spatial width of these averages is plotted in Fig. 8 B. 
As discussed above, spatial width, or FWHM, was well 
defined in the average spark data, independendy of the 
direction in the x—y plane, and there was no significant 
difference between the average values of FWHM in two 
orthogonal directions (x and y dimension). This sym- 
metry is inconsistent with the possibility that a spark is 
originated by an extensive source of size greater than 



the limit of resolution of the microscope. It is not incon- 
sistent with an asymmetric source, constituted by multi- 
ple channels, provided that its size remains below the 
limit of spatial resolution. 

Another intriguing property of these events is re- 
vealed by the study of widths. As the plots in Fig. 8 B 
show, there is no obvious difference in the evolution of 
width among sparks of different rise time, and the plots 
show no visible inflection at the times of the peaks. If there 
is no difference in width among the groups, it follows 
that release termination has no consequences for spark 
width. Models of spark production predict some in- 
crease in width associated with release termination, 
superimposed on the increase mandated by simple 
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Figure 8. Properties of sparks of different rise times. (A) Ampli- 
tudes (i.e., maxima of F/F 0 ) at successive values of t (time from 
first detection) for averages of sparks of different T (coded as in- 
dicated in B) . Note that the early course of fluorescence is similar 
for all averages (i.e., independent of 7). (B) FWHM measured at 
different times on the same averages. The values are similar and 
evolve similarly with ( for all averages. (C) Dependence of spark 
parameters with T. Peak fluorescence (green trace) grows with 
T. Peak flux (red) does not have a clear dependence on T. Flux 
at time of peak fluorescence (circles) decays monotonically as 
T increases. 
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diffusion (a phenomenon named "post-peak expan- 
sion" in Zhou et al., 2003). It can also be seen that the 
increase in width is essentially a linear function of time, 
in every case. Simple diffusion instead predicts an in- 
crease proportional to the square root of time. A possi- 
ble reason for these discrepancies is that the theory in 
Zhou et al. (2003) assumes complete termination of 
Ca 2+ release at the peak and a constant current during the 
postulated release time. As discussed below, neither con- 
dition is likely to be satisfied by the events observed here. 

Peak fluorescence and flux in these averages of sparks 
of different rise time are summarized in Fig. 8 C. Note 
first that spark amplitude (i.e., normalized increase 
in fluorescence at the time of the peak, in green trace) 
increases with T. This property was already presented 
in Fig. 8 A. In its regard, we noted that the early, rapid 
rise in fluorescence was similar for all groups of sparks, 
regardless of their T, implying that release flux was simi- 
lar, and maximal, for all averages at this early time. Peak 
flux, plotted in red, is indeed similar for all groups, and 
so is peak current (not depicted) . 

In Fig. 7 B (illustrating the group of sparks with a 
common T = 11.2 ms) , it can be seen that the flux 
peaked at 5.6 ms, before the time of peak fluorescence. 
The peak of flux was reached at 5.6 ms in all other cases, 
implying that flux decayed during the rising phase of 
the spark. This is demonstrated in Fig. 8 C, where circles 
plot flux at the time of peak Tfor all five average sparks. 
That flux at the time of peak F is lowest for the sparks of 
longest rise time suggests that flux decreases monotoni- 
cally during the rising phase. 

These observations regarding the evolution of [Ca 2+ ] cyto 
and flux have important mechanistic implications. 
The main observations are as follows: (a) F(i), hence 
[Ca 2+ ] cyto (£), dF/dt, and Flux(£) are approximately the 
same for all groups, up until their respective peak 
times ( Tj) . (b) After an early maximum, which occurs at 
5.6-11.2 ms, flux decays substantially during the ris- 
ing phase of the sparks. Observation a suggests that 
[Ca 2+ ]sR(0 is also very similar for these groups, up until 
their respective peaks. 

These data, plus observations in previous studies, pro- 
vide constraints for interpreting the underlying chan- 
nel behavior. Ca 2+ release current is equal to NxPoxi, 
where N*Po represents the number of open channels in 
the couplon, and i represents the unitary channel cur- 
rent. Flux and current are proportional; therefore, a 
decreasing flux or current during the times Tj requires 
the decay of i, N*Po, or both. 

Previous works provide abundant evidence that i de- 
cays (because local depletion, inferred from simulations 
and demonstrated by blinks, implies a decrease in the 
Ca 2+ gradient) . NxPo is also expected to decay, as SR de- 
pletion has a well-established closing effect on channels 
(Sobie and Lederer, 2012). A satisfactory quantitative 
description of sparks and blinks in terms of the above 



factors, however, does not exist (the article by Sobie and 
Lederer, 2012, discusses the difficulties that any attempt 
at such description will find) . The present observations, 
showing that sparks of very different rise times have 
essentially indistinguishable early stages, set additional 
constraints and challenges. At least in an approximate 
and average sense, it can be affirmed that both the num- 
ber of open channels and the unitary current are evolv- 
ing similarly (at times before Tj) in the sparks studied 
here. If that is the case, a justification for the different 
rise times is not evident. If rise times are associated with 
the time when most channels are open (as assumed in 
models of Shen et al., 2004, and Wang et al., 2004), the 
observation would imply that the open times are not 
controlled by [Ca 2+ ]sR as is widely believed, or that the 
control is exerted, but different groups of channels have 
different sensitivity to it. The latter interpretation is sup- 
ported by indications that the [Ca 2+ ]sR level at which 
sparks terminate is variable among couplons, but is 
highly constant for a given individual release site (Zima 
et al., 2008a). 

Other simulations suggest that the duration of the 
rising phase of sparks is determined by the interplay of 
depletion and the varying number of open channels 
(e.g., Sobie etal, 2002). Certainly, different rates of deple- 
tion, associated perhaps with different volumes of SR 
cisternae or different refilling rates (Picht et al., 2011), 
could result in fluorescence peaks at different times. 
The difficulty with such scheme is that it would predict 
a gradual divergence of F between the groups with dif- 
ferent rates of depletion, rather than the sharp depar- 
ture at Tj observed here. 

In conclusion, the present observations add to the 
current picture of control the possibility of other, less 
deterministic mechanisms. One is closure of groups of 
channels triggered by depletion at very different [Ca 2+ ]sR 
levels. A second possibility is that control by [Ca 2+ ]sR is 
weakened and consequendy sparks terminate for rea- 
sons other than local SR depletion (a conclusion consis- 
tent with the observation of long-lasting sparks of varied 
durations but constant [Ca 2+ ]sR; Zima et al., 2008a). The 
latter could be a peculiarity of atrial cells, consistent 
with other indications that the feedback mechanisms 
that control Ca 2+ release are weakened in these cells. 

In-focus sparks have greater correlations between 
morphometric variables 

As demonstrated by 2-D histograms in Fig. 3, the mor- 
phometric parameters a, FWHM, and Twere positively 
but poorly correlated for sparks detected in 3-D (left- 
side panels) . When the analysis was restricted to sparks 
in focus, the correlations increased in every case. 

A positive correlation between rise time and ampli- 
tude is in contrast with previous observations. Indeed, 
in cardiac myocytes, Wang et al. (2004) noted a quantized 
distribution of rates of rise of fluorescence, which they 
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attributed to discrete numbers of participating channels 
(n q ) and a reduction of rise time for sparks of high n q , 
which was interpreted as evidence that large currents 
reduced spark duration by closing channels. In cardiac 
myocytes, the correlation between amplitude, rise time, 
and width was weak (albeit positive; Shen et al., 2004). 
In skeletal muscle of amphibians, an absence of corre- 
lation between these variables was noted (Klein et al., 
1999; Lacampagne et al., 2000). Rfos et al. (1999, 2008) 
found a small and not significant negative correlation, 
which upon a more refined analysis yielded significant 
correlation in two distinct regions: a region of low T where 
the correlation was positive, followed by one of signifi- 
cant negative correlation at longer T. Both an absence 
of correlation and a negative one require some sort of 
negative feedback that turns off channels when local cy- 
tosolic [Ca 2+ ] increases (or [Ca 2+ ] in SR lumen decreases). 
In the present case, however, a positive correlation was 
found. That it became more significant when the analysis 
was restricted to sparks in focus increases confidence in 
the finding. The positive correlation suggests that the 
negative feedback mechanisms on channel openness, 
implied by the aforementioned work, may be less impor- 
tant in atrial cells or under the conditions of the pres- 
ent experiments. The large flux and release current 
calculated for the observed sparks is also an indication 
of increased propensity of channels to open, as is the 
high frequency of these events. 

The positive correlation observed between rise time 
and width is as predicted in spark models (for example, 
Jiang et al., 1999; Zhou et al., 2003). Finally, the positive 
correlation between amplitude and width (which follows 
mathematically when the other two correlations are 
positive) can be understood as the consequence of both 
variables increasing with increasing rise time (depen- 
dencies that were both found, as documented in Fig. 8) . 
These simple positive correlations are predicted by 
models in which width and amplitude are determined 
largely by the duration of the underlying Ca 2+ release 
current (which is positively correlated with rise time) , 
rather than those in which longer durations of release 
are associated with smaller clusters of channels or lower 
unitary current. The latter associations are consequences 
of negative feedback mechanisms affecting channel open- 
ness. Therefore, these correlations are again indications 
of weakened negative feedback in the present case. 

These cells were subjected to membrane permeabi- 
lization, which may have resulted in the increased 
spark size and frequency, in a way similar to the in- 
crease of spark size and frequency observed in cut ver- 
sus intact skeletal muscle fibers (Baylor et al., 2002; 
Chandler et al., 2003) and in chemically permeabilized 
versus intact atrial myocytes (Sheehan et al., 2006), or 
the appearance of spark-like events in mammalian skel- 
etal muscle cells when they are peeled or saponized 
(Kirsch etal., 2001). 



On the other hand, this weakened negative feedback 
could be a feature of atrial myocytes, critical for Ca 2+ 
release during excitation-contraction coupling (ecc) . 
The lack of t-tubules in atrial myocytes requires a dis- 
tinctly different mechanism of Ca 2+ release during ecc. 
Although in ventricular myoctes all RyR clusters are 
found in close physical association with a dihydropyr- 
idine receptor (DHPR) Ca 2+ channel in the surface 
membrane because of the extensive t-tubular system, in 
atrial myocytes only the RyRs of the junctional SR in the 
cell periphery are organized in a similar arrangement. 
Consequently, action potential-induced opening of 
DHPRs provides the required trigger Ca 2+ for CICR si- 
multaneously and homogeneously throughout an en- 
tire ventricular myocyte, but not in atrial cells. In atrial 
myocytes, DHPR opening triggers release initially only 
in the cell periphery, which leads to an elevation of sub- 
sarcolemmal cytosolic [Ca 2+ ] that initiates regenerative 
and propagating CICR from nj-SR. This centripetal 
propagation of CICR from nj-SR is reminiscent of car- 
diac Ca 2+ wave propagation (Berlin, 1995; Huser et al., 
1996; Kockskamper et al., 2001; Mackenzie et al., 2001; 
Woo et al., 2002; Sheehan and Blatter, 2003). Ultra- 
structural studies (Kockskamper et al., 2001) showed 
that nj-SR membranes do not form the narrow diadic 
clefts that would allow [Ca 2+ ] to rise rapidly on the cyto- 
solic face of the RyR. Thus, to achieve reliable and ro- 
bust Ca 2+ release throughout the entire atrial myocyte, 
other mechanisms need to be postulated. As demon- 
strated in the present study, the relative large Ca 2+ flux 
underlying atrial sparks, the calculated large numbers 
of channels involved, the enhanced propensity of RyR 
channels to open, the extended spatial spread (see also 
Blatter et al., 1997), and a diminished influence of nega- 
tive feedback mechanisms in this preparation compared 
with ventricular myocytes all contribute to a situation 
where lower [Ca 2+ ] levels are required to trigger CICR 
and propagation of activation during ecc. 

Conclusions 

We have demonstrated a method to image sparks that 
are in focus, based on scanning fluorescence in three 
spatial dimensions. The sparks found to be in focus had 
substantially greater amplitudes than the population 
of all detected sparks and, at variance with these, 
distributed with a broad mode (at ~0.7 F 0 ). Properties 
of sparks in focus were determined both with accuracy, 
because the out-of-focus error was removed, and preci- 
sion, given the large number of events recorded. Among 
the main morphometric parameters, amplitude was 
substantially different in the sparks in focus, whereas 
average width and rise time were not. The Ca 2+ release 
flux and current underlying sparks in focus were large, 
requiring sources of 20-30 channels. The flux peaked 
early and then decayed during the rising phase of the 
spark, which is consistent with existing evidence of 
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Figure A1 . The three-point criterion for sparks in focus. (A) The spa- 
tial (z) profile of fluorescence of a spark, represented as a Gaussian 
function, 
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of cr = 1.4 um, centered at £ = 0 (black), 0.825 (red), or 1.65 (green). 
Segments ai-a 3 , bi-b 3 , and Ci-c 3 join the intersections of these curves 
with die lines at Zi ( — 1 um) and Zs (+1 um). The difference between 
the ordinates at these intersections increases as die central abscissa f 
moves away from z = 0. (B) S, absolute value of the difference between 
fluorescence at z = Zi and z = Z3, normalized by the value at z = Zg (ac- 
cording to Eq. Al), and plotted versus (, the central position (zvalue) 
of the Gaussian spark. Note diat sparks of cr = 1.4 um will have S < 0.2 
(or 20%) when their central zvalue satisfies the inequalities, —0.37 < f 
< 0.37, in other words, when the origin of release is within a spatial slice 
of height 2 x 0.37 (or 2D) .(G) Generalization of the calculation in B 
to include cr as a variable. S, calculated by Eq. Al , is plotted versus f and 
FWHM = (2 x ln2) 1/2 <r. This function works as a threshold criterion 
(hence it is labeled C3) as follows: given a set thickness 2D of the slice 
of space deemed to be in focus, an experimentally recorded spark 
will be in focus if its calculated 5 is <C3 at the set value of £)and the 
FWHM measured for the individual spark. In the graph, the line 
plots are intersections of the C3 (D, FWHM) surface by planes of 
constant C3 (black), D (red), or FWHM (green). 



significant local depletion of the SR, but could also be 
helped by progressive closure of channels in the cluster. 
The magnitude and evolution of flux were similar in 
sparks of different rise time. This feature indicates that 
the termination of Ca 2+ release was not robustly con- 
trolled by variables associated with flux, including free 
cytosolic and free SR calcium concentrations, either be- 
cause different channels responded to these variables in 
different ways, were controlled by variables other than 
[Ca 2+ ] cyto and [Ca 2+ ]sR, or closed at random. The large 
sparks, large fluxes, and calculated large numbers of 
channels involved define a state of enhanced proclivity 
of channels to open, reflecting a diminished influence 
of negative feedback mechanisms in this preparation. 
The lack of association between flux and spark rise time 
can also be taken as an indication of a general deficit in 
negative controls, underlying the abundance of local 
events in these atrial cells. Because atrial myocytes lack 
a t-tubular system and have to rely on robust cell-wide 
propagating CICR for the activation of the contractile 
machinery, this deficit in negative control constitutes 
an advantage for ecc in these cells. Other properties of 
local events in focus, including spatial associations, event 
propagation, and location-specific aspects of Ca 2+ release, 
can be studied advantageously with 4-D scanning. 

APPENDIX 



A robust procedure for identifying in-focus sparks 

The purpose of this appendix is to describe and illustrate 
a procedure to decide, based on comparison of images at 
three values of the z coordinate, whether a spark is in 
focus. The procedure and associated criterion are called 
"three-point" because they use asymmetrically the mea- 
sures of amplitude at the different z positions. In contrast, 
the simple criterion, whereby a spark is deemed in focus if 
its amplitude at z% is greater than at the other two z posi- 
tions, uses the measures at Z\ and z 3 symmetrically, thus not 
taking advantage of all the information available. 

The procedure is illustrated in Fig. Al. A spark is repre- 
sented by a Gaussian function of the vertical coordinate z. 
The function has a standard error (cr) of 1.4 um, which 
makes it representative of sparks in the present database. 
We will refer to this function as /i(z,1.4). Gaussian sparks 
with identical cr are depicted centered at two other posi- 
tions on the z axis, separated by 0.825 um — an arbitrary 
distance of convenience. Three vertical lines are traced at 
positions Zj ( — 1 pm) , zg (0 pm) , and Z3 (+1 pm) , represent- 
ing the three planes of x-y-z scanning. The segments ai-a3, 
bi-b3, and C1-C3 are the intersections of the functions with 
the planes z\ and z 3 . Their vertical spans are therefore the 
differences between the spark intensity measured at the 
two extreme values of z. It can be seen that this difference, 
as a fraction of the central (z^ = 0 pm) value, increases 
steeply as the spark moves away from the central plane. 
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The absolute value of the difference, as a fraction of 
the value at zg, can be calculated as: 

g = ft(-g-L1.4) -&(-£+ 1,1.4) 
AK.1.4) 

where £ is the separation between the focal plane and 
the center of the displaced spark. 8 is represented as a 
function of £ in Fig. Al B. For absolute values of f <0.37 pm, 
5 is <0.2. In other words, when the difference between 
amplitudes at the two extreme z planes is <20% of the 
central amplitude, events of cr = 1.4 pm should be re- 
stricted to a slice 0.74-pm thick, centered by the focal 
plane at z = 0 pm. 

Inspection of Eq. Al and Fig. Al shows that 8 decays 
monotonically as £ decreases, a property that was used 
to advantage to generate a scalable criterion for identi- 
fying in-focus sparks. Let 8, calculated by Eq. Al for a 
set of values of £, be less than that calculated for a refer- 
ence value of £, say, D. Because of the monotonic rela- 
tionship, it follows that all f that satisfy the inequality 
8 (£) < 8 (D) will be less than D. An equivalent state- 
ment is that the origins of sparks at such values of £ 
will be contained within the slice of half thickness D. 
We represent 8 (D), the "criterion" threshold, by the 
symbol C3. The criterion can be made applicable to 
sparks of other spatial widths simply by introducing the 
corresponding value of cr in place of 1 .4 in Eq. Al . C3 is 
therefore a function of D and cr (or the proportional 
quantity (2 x ln2) 1/2 cr = FWHM). C3is represented as a 
function of D and FWHM in Fig. Al C. (Note that C3 
grows monotonically with D, as shown by the projection 
in green, but has a modal dependence on FWHM, shown 
by the projection in red.) Indeed, for very narrow 
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Figure A2. Convergence of the three-point criterion for iden- 
tification of in-focus sparks. Histograms of amplitudes of sparks 
selected as in focus by the three-point criterion, at progressively 
lower values of D (half-depth of in-focus slice). Values of D are 
indicated on the graph. Note that the features of the histograms 
vary little at D < 0.1 pm. In this work, a spark was considered to be 
in focus if its S, calculated by Eq. Al, fell below the C3 value at the 
measured FWHM on the intercept of the C3 surface (Fig. Al C) 
and the plane at D= 0.1 pm. 



sparks, the 8 value calculated by Eq. Al is small be- 
cause both terms in the numerator are small. For large 
values of FWHM, both terms become large but similar, 
so 8 becomes small again. 

This criterion proved to be superior for several rea- 
sons, which justify their use in the present work. One 
important advantage is that it satisfied expectations of 
convergence. Specifically, when this criterion was ap- 
plied to the same universe of sparks at progressively nar- 
rower depths of the in-focus slice (set by the variable D) , 
it produced groups of sparks tending to a limiting set of 
properties. This is illustrated in Fig. A2 with histograms 
of amplitudes of the groups of sparks selected by the 
three-point criterion at progressively lower values of D. 
As D decreases, the histograms converge to a distribu- 
tion with greater density at high values of the ampli- 
tude. This analysis cannot be taken to a true limit as the 
spark numbers decrease with D, but the convergence 
seemed satisfactory as far as it could be explored. More- 
over, the loss of density at low values of amplitude and 
gain at high values was as expected if the selection pro- 
gressively eliminated out-of-focus sparks. Specifically, 
the histogram obtained with the criterion applied at a D 
of 0.1 pm (green trace in Fig. A2) appears to have limit 
features, while still including a large number of sparks. 
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